Slow Vibrations in Transport through Molecules 
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We show how one can measure the signal from slow jumps of a single molecule between metastable 
positions using a setup where the molecule is fixed to one lead, and one of the coupling strengths is 
controlled externally. Such a measurement yields information about slow processes deforming the 
molecule in times much longer than the characteristic time scales for the electron transport process. 



One of the key ideas in studies of electron trans- 
port through single molecules is the aim to relate the 
properties of the studied microscopic molecule to the 
current flowing through it. Then measuring this cur- 
rent will yield information about the molecule. There 
are many interesting transport phenomena, known from 
larger structures, e.g., semiconductor quantum dots, that 
have been also observed in moleculesi^. However, per- 
haps a feature most specific to the molecular systems is 
the large signature of the mechanical vibrations on the 
transport properties. Such effects include the electron 
shuttling^ and polaronic effects^, e.g., the vibration- 
assisted electron tunneling effect, observed through the 
side peaks in the differential conductance£&2^i at posi- 
tions corresponding to the vibrational frequencies. An- 
other molecule-specific property can be seen when one 
is able to vary the coupling of the molecule to the leads 
between weak and strong coupling limits^. In this case, 
one can quantitatively characterize the different coupling 
strengths, by fitting the experimentally measured con- 
ductance to a fairly generic model describing transport 
through the closest molecular level(s). Such a model re- 
lies on the fact that the molecule is coupled to the leads 
only from one side, allowing one to tune the other cou- 
pling over a wide range. From this fit, one then ob- 
tains four molecule-specific parameters corresponding to 
the two coupling strengths at given positions, an energy 
scale describing the position of the HOMO /LUMO level 
(whichever is closer) and a length scale describing the 
change of the coupling as a function of the distance. 
These parameters can then be used as a fingerprint of 
that particular molecule. 

The typically considered vibrational effects are char- 
acteristic of weak coupling for the electron hopping be- 
tween the leads and the molecule, in which case the vibra- 
tional frequency scales exceed or are of the same order 
as the coupling strength. In such systems, it is essen- 
tial to consider the fairly fast and low-amplitude vibra- 
tions inside a single parabolic confining potential around 
some long-lived metastable position. However, on a much 
slower scale, the molecule may jump between different 
metastable states corresponding to different conforma- 
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tions or positions. Our aim is to discuss in this paper 
how these jumps may be observed and characterized. 




FIG. 1: Schematic illustration of the considered situation: 
Left: A small molecule is connected to the left lead via a 
linker (this particular molecule is from Ref . ITU) . Due to the 
coupling via a single linker, the molecule-linker-lead system 
has multiple metastable configurations, corresponding to dif- 
ferent positions. The distances d r L and djj to the two leads 
are also indicated. Right: A possible potential profile corre- 
sponding to the position of the molecule, shown in the left 
figure by the coordinate x. Due to the directional character 
of the bonding to the leads, due to the solvent, or due to 
the particular atomic arrangement of the left lead, there may 
be a few metastable configurations indicated by the potential 
minima. The hopping between these configurations depends 
on the heights Ui and Ui of the potential barrier. If Ui U2, 
the probability pi to occupy state 1 is larger than the corre- 
sponding probability P2 for state 2. 

Consider a potential energy curve depicted in Fig. ^ 
The horizontal axis could quantify different molecu- 
lar conformations or average positions. The vibrations 
within a single potential well are governed with a fre- 
quency ujk — y/k/m, where k is the spring constant de- 
scribing the potential and m is the mass of the molecule. 
Brownian motion of the particle within this potential well 
at the temperature T will then result into vibrations with 
amplitude s ~ ^/fc^T/fc. The amplitude of such vibra- 
tions is much smaller than the distance between succes- 
sive potential minima, and hence it is at most of the 
order of a few A. We can use this as an estimate for the 
frequency Uk- At room temperature, for the case of a 
molecule with mass of the order of 1000 . . . 10000 m v , 
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we would thus get u>k > 0.1 ... 1 THz. These vibrations 
are damped by a friction force described by the charac- 
teristic rate 7/ = crjR/m, where c is a constant of order 
unity depending on the shape of the molecule (c = 6n 
for a spherical molecule) , R characterizes the size of the 
molecule, and 77 is the viscosity describing the molecule 
environment. For a spherical molecule of size R ~ lnm, 
mass as given above and using the viscosity of water, 
r/=l g/(ms), we would thus get 7/ ~ 10 THz. Note that 
in practice, the effective viscosity of the solvent depends 
also on the molecule itself and thus this number should be 
used as indicative only. The jumps between the different 
potential minima have a much lower rate than the small- 
scale vibrations. This rate is described by the Arrhcnius 
law* 



7 = v exp(- 



k R T J 



(1) 



where U\ describes the height of the potential barrier (see 
Fig. 0, v = in the overdamped limit jf 3> oj c and 

v = ojfc/(27r) for 3> 7/- Here oj c describes the width 
of the potential barrier, and is of the same order as uik- 
With the above estimates for the frequencies u>k, u c and 
7/, the prefactor v thus ranges from GHz to THz. How- 
ever, the exponential factor makes the jumps between 
different minima much less frequent. Assume for exam- 
ple a potential barrier height of U w 0.5 eV 15 . At room 
temperature, we would then get 7 = z/exp(— 20); ranging 
between Hz and kHz. This is close to the characteristic 
scale in which the measurements on the molecules are 
made and indeed such measurements^ showed large fluc- 
tuations in the measured conductance, clearly connected 
to the presence of the molecule. 

The distance-dependent linear conductance G = GqT 
through a single molecular level can be described by the 
Breit-Wigner formula^, 



T = 



Tr.T 



L L R 



el + {T L + T R y/A- 



(2) 



Here Go = 2e 2 /h, ei is the energy of the closest molecular 
level to the metal Fermi energy (i.e., LUMO or HOMO, 
whichever is closer) assuming it has an appreciable cou- 
pling to the leads, and Ti, and Tr characterize the cou- 
pling to the left and right leads, respectively. The level 
ei may be degenerate - this degeneracy would only tune 
the effective coupling strengths Tl and Tr compared to 
the non-degenerate case. For simplicity, we neglect in- 
teraction effects. This assumption still captures the es- 
sential physics in the strong-coupling regime where the 
coupling energy Tl + Tr exceeds the thermal energy and 
thus describes the lifetime of the level. Moreover, addi- 
tional molecular levels may be considered, but their con- 
tribution shows up mostly to slightly rescale the coupling 
constants^. 

Consider now what happens if the molecule is con- 
nected to one of the leads, say left, thus fixing the av- 
erage r^. Assume furthermore Tl < ei. The average 



coupling to the right lead depends on the distance d,R 
between the molecule and the furthermost atom of this 
lead through Tr = T R exp(—K,d R ), where k depends on 
the solvent and on the molecule/lead materials^. For 
Tr <C ei, decreasing cLr will increase the conductance. 
However, when the right lead is close enough, Tr may 
exceed the level energy e±. In this case, the conductance 
shows a maximum at Tr — \J 4e 2 + T 2 L ~ 2e\ and fur- 
ther decrease of d,R leads to a decrease in the conduc- 
tance. This type of a model was employed to explain the 
observed conductance-distance curve in Ref. (TJ with a 
quantitative agreement between the theory and the mea- 
sured average conductance. 

Consider now the fluctuation of this conductance, due 
to the slow hoppings of the molecule between different 
average positions. Such hopping corresponds to a ran- 
dom telegraph noise in a time-dependent signal. Let us 
denote the average distance between the right lead and 
the molecule by d (average meaning averaging over the 
different positions of the molecule corresponding to the 
given positions of the leads). Let us furthermore choose 
the coupling strengths corresponding to this average po- 
sition to Tl and Tr = T R exp(— nd). Then, the fluctua- 
tions of the position around this average position can be 
characterized by the values {Sd 1 , —aSd 1 }, i indexing the 
different potential minima, and the two numbers corre- 
sponding to the deviations of the distance to the left and 
right leads, respectively. In a typical case, one could ex- 
pect that if the molecule moves further from the left lead 
(Sd > 0), it comes closer to the right lead (as in Fig.Q. 
This would thus correspond to a positive Cj. However, 
for certain situations it may be possible to increase the 
distance to both leads - this would be described with a 
negative Cj. With these deviations, the couplings change 
to = T L exp(K,5d l ) and T R = T R exp(— KCiSd 1 ). Note 
that choosing the same K for both Tl and Tr does not 
mean a lack of generality, as a possible difference in the 
two k's can be included to scale Cj. 

The transmission averaged over the positions of the 
molecule is 



(T)=J2PiTi = 



+ (Ti + ryv4 



(3) 



Here pi is the probability for the molecule to be in the po- 
sition/configuration i. Thus, already the average trans- 
mission (T) depends on the amplitude Sd 1 of fluctuations 
(see Fig. 0). However, such a dependence is difficult to 
see in (T), as a similar behavior could be observed also 
without vibrations, but with a slight rescaling of T^m. 

The variance of the transmission values due to these 
slow fluctuations is var(T) = ((T — (T)) 2 ). Assuming we 
can neglect the electronic noise (see below) which also 
shows up as a temporal variation of the current, this 
var(T) would vanish without the vibrations. In general, 
var(T) depends on d, the separation of the leads. How- 
ever, if we assume that the positions of the metastable 
states with respect to the left lead are independent of 
d, we can separate two limits in the d-dependence. One 
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FIG. 2: Average transmission as a function of coupling Tr for 
different amplitudes of variations Sd. The coupling Tr can be 
varied by varying the distance between the electrodes, and 
the transmission peak details information about the molecule 
and its coupling to the leads. The hopping between differ- 
ent metastable configurations corresponding to different Sd 
affects the transmission: from top to bottom, Sd = (no 
hopping), 5d = 0.5/re and Sd = 1/k. We chose Tl — O.Olei 
and c = 0.5 and describe hopping between two degenerate 
configurations (i.e., with equal probabilities) separated from 
the average position by ±<5d. The dashed lines indicate fits to 
the Breit-Wigner transmission with no account of the fluctu- 
ations, but with an increased Tl and a smaller T* R . 



is the case when d is large, such that Tr -C e%. If also 
Tl *C £i, we can neglect the lifetime of the level (the 
term T L + Tr in the denominator of Tj). Then 



r«2 £1 r^r} 



TlTr , . 
— 2 — exp(/t(l 



c % )bd % ). (4) 



Thus, the transmission probability for each i can be writ- 
ten in a form 7$ — TTj^ where T is independent of the 
random hoppings, but depends on the position d, and 
depends on the hoppings, but not on the position d. 
In this case, we may express the relative variance as 



^2 



= var(T) _ vargFji) 
J2iPi |exp[«5d f (l 



(5) 



Thus, this quantity no longer depends on the exact value 
of d, as long as 1^ + Tr <C 2e%. 

The same happens in the opposite limit, Tr 3> 2ei. In 
this case, we can neglect all other terms but Tr from the 
denominator of the transmission and 



r„> 2ei 4r* £ _ r L 



exp(«;(l + Cl )W l ) = TT: 



/2- 



(6) 



The relative fluctuations <tt again follow Eq. (J5J, with 
the only exception that the sign of each Ci is reversed. 



If the sign of c, is predominantly positive, err m the 
case Tr » 2ei will be larger than in the case T <c 2ei 
and vice versa for a predominantly negative ci. Thus, 
we can sketch the rough behavior of <jt as a function 
of the distance d (assuming c > 0): At first, when the 
leads are far apart, ax stays mostly constant. When 
Tr becomes of the order of 2e±, ot starts to increase 
with d, until saturating into another constant value at 
» 2ei (see Fig. 01. Such a behavior holds as long as 
the metastable positions of the molecule are unaffected 
by the right lead. The latter type of a mechanical effect 
would show up also in the average conductance curves 
if the right lead changes the potential landscape seen by 
the molecule. This was probably observed in Ref.O, but 
only when Tr was already much larger than ei. 
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FIG. 3: Relative variance in the measured conductance 
curves, observed as time-dependent fluctuations. The two 
solid lines are for hopping between two degenerate configura- 
tions with nSd = 0.5 (bottom) and nSd = 1 (top), and the 
dashed curve represents the case with four degenerate config- 
urations, at Sd = ±0.5/k and Sd = ±l//t. We chose c = 0.5 
for each curve. 

To explore this behavior explicitly, let us con- 
sider a simple two-position model with the positions 
{Sd/2, ~c5d/2} and {—5d/2, cSd/2}. In this case, we get 
a fairly simple expression for ctt m the limit Tl <C £i, 



(JT 



4ej sinh( K (l - c)Sd/A) + T R sinh(«(l + c)5d/4) 



—2 



4ej cosh( K (l - c)Sd/A) + T R cosh(K(l + c)Sd/A) 



In the limit Tr <C 2ei, this gives 

aT r «S 2ei | tanh(«(l - c)<M/4)|, 
and in the opposite limit Tr 3> 2e±, 



ax 



|tanh(K(l + c)M/4)|. 



(7) 

(8) 

(9) 



These limits follow the qualitative discussion above. 

Apart from hopping between different positions, in 
some cases one may also envisage the molecule to hop 
between different conformations on the slow time scales. 
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Such a change in the conformation in general may lead 
to a change both in the energy level t\ and in the cou- 
pling strengths Y^m. This behavior can be illustrated 
by considering the simplest case of hopping between two 
conformations corresponding to the energies ei ± <5e/2 
and couplings T L ± ST L /2, T R ± 8T R /2. The relative 
variance gt now depends on the relative magnitude of 
these changes: if 5T l / r /(Tl + T R ) ;§> Se/ei, the behav- 
ior is analogous to that discussed above. In the opposite 
limit of large Se, the relative variance of the conductance 
values g t — var(T)/(T) 2 is given by 



(<5e 2 + (r L + Y R y + 4ef) 2 ' 



(10) 



Thus, the relative variance is largest when the couplings 
are much smaller than the level energies, and it decreases 
as either of the couplings is increased. A similar conclu- 
sion can be drawn for the general case with many differ- 
ent conformations, along the same arguments as above. 

There are a few experimental constraints for the ob- 
servation of the predicted behavior in the fluctuations, 
characterized by the different time scales in the problem. 
An easily satisfied condition is that the measurement 
time r m should exceed the time scales l/u)k, fr/^L/Ri 
r e = e/ (I) characterizing the individual charge transport 
processes (typically between ps and ns) by a few orders 
of magnitude. Here (/) is the average current through 
the molecule. In this limit, shot noise yields a contribu- 
tion ~ e/(r m (/}) = r e /r m to the relative variance and 
can hence be neglected. The same applies for the ther- 
mal noise provided that k R T / {eV)T e /T m <C gt, where 
V is the bias voltage applied over the sample. Another 
natural condition is that the time scale T var for the vari- 
ations made in the structure (like changing the distance 
between the leads) should be longer than r m and the time 
scale Thops = 1/7 for the slow changes in the configura- 
tions. To obtain a relative accuracy p for the measured 



variance, one has to measure at least ~ 1/p 2 points and 
therefore r var /r m > 1/p 2 . 

If there are only a few metastable configurations in the 
problem, and the time scales for hopping between them 
is longer than the measurement time, one may be able to 
measure the information about them already by follow- 
ing the telegraph noise in the average transmission as a 
function of time. However, for many configurations, or if 
at least some of the hopping time scales are smaller than 
T m , it is better to measure the relative variance. When 
r m and Thops are well separated, the measured variance 
in the signal will be proportional to 



var(G) m = var(G), 



min(T m , T hops ) 
max(T m , Thops) 



(11) 



Here var(G) m is the measured variance and var(G) c = 
Gq var(T) is the variance calculated above. In the case 
when there are multiple time scales describing the slow 
fluctuations, and the measurement time is between these 
scales, the measured variance will be independent of r m , 
characteristic for flicker noise. 

Summarizing, in this paper we predict that the differ- 
ent metastable atomic configurations in molecular junc- 
tions have a considerable effect in the measured conduc- 
tance, as the time scale of typical conductance measure- 
ments is of the same order as the time scales for the jumps 
between the different configurations. We utilize a simple 
Breit-Wigner model to illustrate this behavior and show 
that such variations lead to a fairly universal behavior in 
the relative variance of the measured conductance values 
as one of the coupling constants between the molecule 
and the leads is controlled. 
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